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The yield ratios of neutron-proton (R(n/p)) and 3H-*He (R(*H/*He)) with reduced rapidity from 0 to 0.5 
were simulated at 50 MeV/u even-even 2°-°°Ca + "Ca, even-even Î5-78Ni + 58Ni, and 19°-139gp (every third 
isotopes) + '!?Sn for full reduced impact parameters using the isospin-dependent quantum molecular dynamics 
(IQMD) model. The neutron and proton density distributions and root-mean-square radii of the reaction systems 
were obtained using the Skyrme-Hartree-Fock model, which was used for the phase space initialization of the 


projectile and target in IQMD. We defined the unified neutron skin thickness as ARnp = (r2)* ae (y (2, 
which was negative for neutron-deficient nuclei. The unified AR,np values for nuclei with the same relative 
neutron excess from different isotopic chains were nearly equal, except for extreme neutron-rich isotopes, which 
is a type of scaling behavior. The yield ratios of the three isotopic chain-induced reactions, which depended on 
the reduced impact parameter and unified neutron skin thickness, were studied. The results showed that both 
R(n/p) and R(*H/*He) decreased with a reduced impact parameter for extreme neutron-deficient isotopes; 
however, they increased with reduced impact parameters for extreme neutron-rich isotopes, and increased with 
the A Rnp of the projectiles for all reduced impact parameters. In addition, a scaling phenomenon was observed 
between AR, and the yield ratios in peripheral collisions from different isotopic chain projectiles ( except for 
extreme neutron-rich isotopes). Thus, R(n/p) and R(?H/*He) from peripheral collisions were suggested as 
experimental probes for extracting the neutron or proton skin thicknesses of non-extreme neutron-rich nuclei 
from different isotopic chains. 
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I. INTRODUCTION 


of neutron-proton [31-39], and 3H/*He [39, 40], have been 


The exotic nuclei far from the (-stability line have attracted 
widespread attention because of their special properties, such 
as extreme N/Z ratios or isospin, significantly different den- 
sity distributions of neutrons and protons (so-called skin or 
halo structures), small valence nucleon separation energy, and 
increased reaction cross-section. Radioactive ion beam fa- 
cilities have been constructed worldwide to produce exotic 
nuclei and explore their properties. Classical nuclear struc- 
ture models have been continuously developed to comprehen- 
sively study their attractive structural characteristics and im- 
prove the theoretical nuclear potential and equation of state 
(EOS). These include the Skyrme-Hartree-Fock (SHF) [1- 
4], shell model [5-7], relativistic-mean-field [8-11], and so 
on [12, 13]. One of the most fundamental properties of ex- 
otic nuclei is their proton and neutron density distributions. 
The proton (charge) density distribution can be measured with 
high accuracy by the elastic electron-scattering and muonic 
atom X-rays methods [14]. Measurement results of the neu- 
tron radius vary significantly owing to the complex strong 
interaction between nucleons, including the hadron scatter- 
ing [15-17], giant dipole resonance [18, 19], and antiprotonic 
atoms methods [20-22]. Some of the experimental probes 
are unsuitable for unstable nuclei; therefore, studies on nu- 
clei radii primarily focus on stable nuclei or nuclei close 
to {-stability. Some isospin- dependent quantities from nu- 
clear reactions, such as collective flows [23-30], yield ratios 


* Supported by National Natural Science Foundation of China under Grant 
No. 11405025. 
t Corresponding author, ytz0110@ 163.com 


studied to explore the isospin physics such as the symme- 
try energy [41, 42] and nucleon-nucleon cross section [43]. 
We used the neutron skin thickness (defined as the difference 
between the root-mean-square (rms) radii of the neutron and 


proton, i.e., ARnp = (r3) (r2)0/), which was calculated 
using the SHF method, and explored the yield ratios of neu- 
tron/proton and 3H/*He within projectile-like rapidities from 
different Ca isotopic projectiles (the collisions were simu- 
lated in the framework of isospin-dependent quantum molec- 
ular dynamics (IQMD) model using nucleon density distribu- 
tions from SHF above for phase initialization of Ca projec- 
tiles) dependent on the neutron skin thicknesses of Ca pro- 
jectiles, which show linearly increasing correlations between 
them [44]. Similar correlations were studied for the proton 
skin thicknesses of neutron-deficient Ni isotopes that exhib- 
ited exponentially declining relationships [45], where the pro- 


ton skin thickness was defined as AR,, = (ry? — (r31, 


In this work, we used the unified neutron skin thickness 
(ARnp = (r3? — (r2)1/?) to combine the neutron skin 
thickness and the proton skin thickness, such that the ARnp 
of a neutron-deficient nucleus was negative in that case. The 
unified neutron skin thickness dependence of the yield ratios 
of the neutron-proton and 3H-?He was studied for Ca, Ni, 
and Sn isotopic projectiles, and the reduced impact parameter 


dependence for the above correlation was also explored. 


Il. METHOD 


The neutron and proton density distributions and nuclear 
radii of the projectiles and targets were extracted using the 
SHF theory with the so-called SkM* parameters [46]. These 


were then used in the IQMD model to generate the phase 
space of the projectile and target. The proper phase space 
must be correctly selected according to the experimental aver- 
age binding energy and rms radius. The time-evolution stabil- 
ity of the phase space should be carefully examined to avoid 
large fluctuations in the average binding energy and rms ra- 
dius. The selected phase space should have a small fluctuation 
in the rms radius, particularly in the first several dozen fm/c 
and for nuclei far from the -stability line. Using these care- 
fully selected initialization phase spaces of the projectile and 
target as input files, further time evolution of the reaction was 
performed within the framework of the IQMD model, which 
is based on the general QMD model [47]. 

The nuclear mean-field potential used in the IQMD model 
can be parameterized as follows: 


U(p,tz) =a (4) +a(4) +30- 


(Pn = Pp) 
Po 


() 


+Cyym ppU TE 


where po is the normal nuclear matter density (0.16 fm~°); 
Ps Pp» Pn are the total, proton, and neutron densities, respec- 
tively; and 7, is the zth constituent of the isospin degree of 
freedom, which is equal to 1 or —1 for neutrons or protons, 
respectively. The parameters a, 3, and y represent the differ- 
ent nuclear equations of state. Cyym is the symmetry energy 
strength owing to the asymmetry of neutrons and protons in 
the nucleus. Cyym = 32 MeV was used in this study. V, and 
Uyak denote the Coulomb and Yukawa (surface) potentials, 
respectively. We used aœ = -356 MeV, 6 = 303 MeV, and y 
= 1.17, which correspond to the so-called soft EOS with an 
incompressibility of K = 200 MeV. Although a soft EOS was 
used without momentum-dependent interactions and yield of 
fragments was affected, the yield ratio remained relatively un- 
changed. 

In this study, the reduced impact parameters of collisions 
of even-even 36-56Ca + 4°Ca, even-even *8-78Ni + 58Ni, 
and 100-139Sn (every third isotope) + !!?Sn at 50 MeV/u 
were simulated using the IQMD model with the soft EOS 
parameterization. The reduced impact parameter is defined 
as bre = b/bmax, Where bmax is the sum of the radii of the 
target and projectile. The yields of light fragments were re- 
stricted within projectile-like midrapidity; therefore, the re- 
duced rapidity of the fragments (y = (y/Yp),,,» Where yp is 
the projectile rapidity) was between 0 and 0.5 to ensure that 
the fragments originated from the projectile. The yield ratios 
of neutron-proton (R(n/‘p)) and 7H—*He (R(?H/*He)) were 
calculated using the yields of neutrons, protons, 9H, and He. 
The yields of the light particles were accumulated from 160 
to 200 fm/c to improve statistics. 


Ill. RESULTS AND DISCUSSION 
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Fig. 1. (Color online) Normalized neutron (solid lines) and proton 
(dashed lines) density distributions of 364%56Ca, 48:5878Ni, and 
100,112,139 calculated using SkM* parametrization of the SHF 
method. 
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Fig. 2. (Color online) Left panel: Mass number dependence of the 
rms radii of neutrons (line + square) and protons (line + circle) (left 
panel) of even-even 36-56 Cg even-even 48—78Ni and 100-139¢h 
(every third isotope) calculated using SkM* parametrization of the 
SHF method. The solid and the dotted lines are the referenced re- 
sults for neutrons and protons, respectively, in which even-even Ca 
isotopes are taken from Ref. [48] calculated with relativistic contin- 
uum Hartree-Bogoliubov (RCHB) theory. The even-even Ni and Sn 
isotopes are extracted from Ref. [49] in the framework of the de- 
formed self-consistent mean-field Skyrme HF+BCS method. Right 
panel: Same as left panel but for the unified neutron skin thickness. 
The upward triangles are our SHF results; the solid lines are the cal- 
culated results based on the referenced neutron and proton rms radii 
of the left panel; and the downward triangles are the experimental 
values extracted from Ref. [50]. 


The normalized neutron and proton density distributions of 
typical neutron-deficient, stable, and neutron-rich isotopes of 


Ca, Ni, and Sn obtained using the SHF method with the SkM* 
parameters are shown in Fig. 1. The proton skin and neu- 
tron skin structures are evident in the neutron-deficient and 
neutron-rich nuclei of medium-mass Ca and Ni, respectively, 
whereas the proton skin structure is difficult to form in the 
relatively heavier element Sn owing to the Coulomb poten- 
tial, as we have already realized. The nuclear radii of the iso- 
topes of the three elements (Ca, Ni, and Sn) related to their 
nuclear mass numbers are shown in the left panel of Fig. 2. 
The squares and circles represent the rms radii of neutrons 
and protons, respectively. The solid and dotted lines indi- 
cate the results for even-even Ca with the RCHB method [48] 
and Ni and Sn with the deformed self-consistent mean-field 
Skyrme HF+BCS method [49], which are close to our results 
for Ni and Sn but deviate significantly for the proton density 
of Ca. The calculated unified neutron skin thicknesses of Ca, 
Ni, and Sn, based on the simulated neutron and proton radii, 
are shown in the right panel of Fig. 2. The upward trian- 
gles and solid lines are our SHF results and the theoretical 
results from the aforementioned studies, respectively, while 
the downward triangles are experimental values of these three 
nuclei extracted from Ref. [50] obtained using the antiproton 
atoms method. Our simulated A Rpp closely matched those of 
other studies, except for extremely neutron-rich Ca, and were 
within the error bars of the experimental values. The exper- 
imental data for the Sn isotopes measured using other meth- 
ods were also compared with the model results in Ref. [49] 
and were consistent with experimental results within the er- 
ror ranges. 

Reference [50] demonstrated the existence of an approxi- 
mate linear correlation between the experimental value of the 
neutron skin thickness of 26 stable nuclei and their relative 
neutron excess I (I = (N — Z)/A). We present the same 
correlation for the Ca, Ni, and Sn isotopes in Fig. 3. The 
simulated unified neutron skin thickness with the SHF theory 
increases linearly with the relative neutron excess, except for 
those of extreme neutron-rich nuclei, and the curves for the 
three isotopic chains almost overlap, except for those of ex- 
treme neutron-rich isotopes. This suggests the possibility of 
scaling behavior for different isotope chains for the correla- 
tion. This interesting phenomenon can further elucidate the 
EOS of asymmetric nuclear matter [51], and more accurate 
theoretical calculations and experimental data on the unified 
neutron skin thickness of isotope chains are required. 


The extracted projectile-like yield ratios R(n/p) (left pan- 
els) and R(?H/*He) (right panels) from the Ca (upper pan- 
els), Ni (middle panels), and Sn (bottom panels) isotopic pro- 
jectiles within different reduced impact parameters are shown 
in Fig. 4. Approximately half of the isotopes are shown in the 
figure to illustrate the changing trend more clearly. The val- 
ues of R(n/p) and R(?H/*He) from different bye are differ- 
ent. They decrease with increasing bre for extreme proton-rich 
isotopes, and gradually change to exhibit an increasing trend 
with the bie with the growing constituent neutron-proton ra- 
tio (N/Z) of the projectile isotopes. Furthermore, larger val- 
ues at the same b,. correspond to projectiles with larger N/Z. 


np 
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Fig. 3. (Color online) Relative neutron excess dependence of the 
unified neutron skin thickness simulated for the Ca (squares), Ni 
(circles), and Sn (triangles) isotopes. 


The fragments at 0 < y < 0.5 are emitted from the projec- 
tile in the overlap zone of the collision, and their constituent 
neutron-proton ratios should be close to that of the projec- 
tile. However, as the proton or neutron skin structure plays an 
increasingly important role with increasing collision param- 
eters, the N/Z of the projectile-like fragment and the yield 
ratios of R(n/p) and R(?H/*He) vary accordingly. Thus, 
the skin structure can affect the yield ratios, and therefore, 
the unified neutron skin thickness dependence of the yield ra- 
tios from different reduced impact parameters for Ca (upper 
panels), Ni (middle panels), and Sn (lower panels) are ex- 
amined in Fig. 5. The yield ratios R(n/p) and R(?H/*He) 
from all the reduced impact parameters both increase with the 
unified neutron skin thickness A Rpp in a monotonic but non- 
linear manner, and the increase is faster for larger reduced im- 
pact parameters. The corresponding nuclei with equal yield 
ratios can also be easily obtained from different reduced im- 
pact parameters. The left panel of Fig. 5 indicates that 4°Ca, 
60Ni, and '!8Sn have equal yield ratios R(n/p) for different 
reduced impact parameters, and they are near the 8 stability 
line. As shown in the right-hand panel of Fig. 5, °°Ni and 
106Sn have equal R(?H/*He), which are neutron-rich nuclei 
with smaller N/Z than those obtained from R(n/p). The Ca 
isotope with equal R(?H/*He) from different reduced impact 
parameters is beyond the Ca scope we studied; however, it is 
definitely neutron-rich. 


Figure 6 and Fig. 7 compare the isotopic chains of the 
projectile-like yield ratios R(n/p) and R(?H/*He) from dif- 
ferent reduced impact parameters that depend on the unified 
neutron skin thickness AR,,, of the projectile. Differences 
are observed in the yield ratios among the three isotopic sys- 
tems in central and semi-peripheral collisions, especially for 
Sn projectiles. However, in peripheral collisions, the yield ra- 
tios from the three isotopic systems almost overlap with each 
other, except for several extreme neutron-rich Ca isotopes. 
This may be due to the inaccurate neutron skin thickness of 
extreme neutron-rich light Ca isotopes calculated using the 
same set of theoretical parameters. Thus, a scaling behavior 
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neutron-rich isotopes) between the projectile-like light frag- AR, (fm) 


ment yield ratios from the peripheral reactions and the unified 
neutron skin thickness of the projectile. In other words, nu- 
clei with the same unified neutron skin thickness would pro- Salon 
duce the same R(n/p) and R(?H/*He) from peripheral colli- A He), 
sions. The scaling behavior can be used to extract the unified 
neutron skin thickness (including neutron-rich or neutron- 
deficient isotopes) of an isotope by experimentally observing 
the mid-rapidity yield ratios of R(n/p) or R(?H/?He) in the 


Fig. 7. (Color online) Similar to Fig. 6 but for the yield ratio 


peripheral reactions from the isotopes of other elements. The IV. SUMMARY 
enveloping curves for R(n/p) and R(?H/?He) for the three 
isotopic chains from peripheral reactions are plotted in Fig. 6 The projectile-like yield ratios R(n/p) and R(?H/*He) 


and Fig. 7 to show the estimated range of skin sizes extracted were studied for even-even 36-56Ca + 4°Ca, even-even 


48—78Ni + °SNi, and !°°-139Sn (every third isotope) + 112Sn 
at 50 MeV/u simulated using the IQMD model with the sam- 
pled initial nucleon densities by applying the SHF method. 
The yield ratios decreased with the increasing reduced impact 
parameter for extreme neutron-deficient isotopes, and then 
gradually increased with the reduced impact parameter for 
extreme neutron-rich isotopes. The yield ratios from different 
reduced impact parameters for isotropic projectiles monoton- 
ically increased with the unified neutron skin thickness A Rpp 
of the projectile. Scaling behaviors may exist between the 
unified A Rpp and the relative neutron excess for different iso- 
topic chains, and between the unified AR,p of the projectile 
and the projectile-like light fragments’ yield ratios from pe- 
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